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Ph! Abstract 

r-| . A measurement of the hadron energy spectrum in inclusive semileptonic B decays is proposed 
^ ■ as a viable method for extracting \Vub\- Compared to the traditional energy spectrum of the 
charged lepton, the hadron energy spectrum exhibits kinematical advantages such as a wider 
^ • energy window and a larger signal branching fraction. It is emphasized that the hadron energy 
. i spectrum method is most suited for symmetric B factories, such as CLEO. The hadron energy 
distribution is calculated in the approach of the Altarelli et al. model and of the heavy-quark 
effective field theory. In both methods, perturbative QCD corrections, the Fermi motion of the 
6-quark in the 5-meson, and the recoil momentum of the S-meson (stemming from the T{4S) 
resonance) are taken into account. We have found excellent agreement between the spectra 
calculated in both methods, especially in the relevant kinematical region below the charmed 
meson threshold. The theoretical error to \ Vub\, which is dominated by the uncertainty of the 
6-quark mass, is estimated to be at the ±12% level. 
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1 Introduction 



A precise determination of the Cabibbo-Kobayashi-Maskawa matrix element V^b is an important 
step for constraining the unitarity triangle. Therefore, it poses a challenge for both theory 
and experiment Traditionally, Vub has been extracted from the energy spectrum of the 
charged lepton in inclusive semileptonic B decays B —>■ Xud-v above charm threshold, i.e., 
for lepton energies E^^ above (m| — mj^) / {2mB) ~ 2.3 GeV [0, 0] 0. As the kinematical 
endpoint for semileptonic B decays is at Eg = {m% — m\) / {2mB) ~ 2.6 GeV, the pure b — *• 
u^~V transition extends over a relatively narrow window of about 300 MeV, accepting only 
a fraction of approximately 10 % of the total sample of charmless B decays. Theoretically, 
the lepton energy spectrum may be calculated from first-principle QCD. Using the tools of the 
operator product expansion (OPE) and heavy-quark effective field theory 0, one can construct 
a systematic expansion of the lepton spectrum in powers of A/m^, where A is a typical low- 
energy scale of QCD 0. However, the relevant lepton energy window from which \Vub\ can 
be extracted lies to a large extent in the so-called endpoint region which extends from the 
parton model maximum at Ei = mb/2 up to the hadronic maximum at {m\ — m^)/(2mB). It 
is therefore described by genuinely non-perturbative contributions. The difficulty arises from 
the fact that, close to the partonic endpoint, the expansion parameter is no longer A/rrib, 
but A/{2mb — Eg). Thus, the theoretical prediction becomes singular when the lepton energy 
approaches the parton model endpoint; formally, these singularities manifest themselves in 
delta functions and derivatives of delta functions concentrated at the partonic lepton energy 
endpoint mb/2. In addition, this region is plagued by large perturbative Sudakov-like double 



logarithms |10|, ^ as well as by small instanton effects |T2|. Therefore, in a large part of the 
lepton energy region in which the extraction of Vub is kinematically possible, a full resummation 
of Sudakov-like double logarithms and power corrections ||13[ becomes necessary. In view 
of the theoretical and experimental difficulties just mentioned, it is desirable to look for other 
methods. 

As an alternative, we propose to extract Vub from the hadron energy spectrum in the inclusive 
charmless semileptonic B decays B — »• Xu^v- Simple kinematical considerations support why 
this proposal is viable. As the charmed final state threshold is at the D-meson mass, the 
maximal hadron energy window for charmless semileptonic B decays is given by the range 
'^TT < Ehad < ^D'l this window is much wider than the corresponding kinematical window 
for the lepton energy distribution discussed above, and as we will see later, a much larger 
fraction of the B — >• X„£z/ events becomes accessible, leading to improved statistics. Of course, 
the theoretical problems addressed in the discussion of the lepton endpoint spectrum are also 



present in principle in the hadron energy spectrum [|I4|, |T5[, but to a much lesser extent in 
the kinematical region relevant for the extraction of Vub- Indeed, up to perturbative QCD 
corrections, the region around Ehad = mb/2 is also fully dominated by non-perturbative effects, 
quite in analogy to the lepton energy endpoint region. Fortunately, this region is well above 
the charmed hadron final state threshold lying outside the region we are interested in. At the 
lower end of the hadron energy, Ehad ~ 0, the hadronic mass ranges over a narrow window 
< mj^^^ < E'^^^, hence, the OPE breaks down again. This region can be avoided by applying 
a lower cut to the hadron energy; we choose the value 1 GeV. This relatively high lower-cut has 



^ We will discuss the impact on the determination of \Vub\ from the recent CLEO measurements of the 
corresponding exclusive channels B (tt, p, lu) £i/ 0, |^ later. 
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the advantage that a wide range of invariant hadronic masses contributes to the hadron energy 
spectrum; quark-hadron duahty, which we imphcitly assume in our treatment for the inclusive 
B Xuiu decay, is then expected to work well. Even after this cut at 1 GeV is made, we are 
still left with an ample hadron energy window ranging from 1 GeV to rriD. These features in 
principle make investigations with hadron energy spectra more reliable than with lepton energy 
spectra. 

Experimentally, the hadron energy spectrum in semileptonic B decays may be measured 
schematically as follows: Working at the T{4S) resonance, which decays into BB, one requires 
one of the 5-mesons to decay semileptonically and the other one hadronically. In the case of a 
symmetric 5-factory, like CLEO, the energy of the hadrons stemming from the semileptonically 
decaying 5-meson can be obtained by measuring the total energy of all the hadrons in the final 
state and then subtracting mx(45)/2. In case of asymmetric B factories, the hadron energy 
spectrum is harder to measure. One way is to reconstruct in a first step the whole T(4S') 
decay in its rest frame and then perform the analysis just described for the symmetric case. 
To perform the corresponding boost, one has to measure precisely both the energy and the 
momentum of each final state hadron, which requires accurate particle identification. 

The remainder of this paper is organized as follows. In section 2 the calculation of the hadron 
energy spectrum is presented. We utilize two methods to account for the bound-state effects. 
The first one, discussed in section 2.1, uses the Altarelli et al. (ACCMM) model [0 of the 
5-meson and the second one, presented in 2.2, is based on the heavy-quark effective field theory 
(HQET) [Q. In both treatments, the perturbative 0{as) corrections are taken into account. In 
section 3 the extraction of \Vub\ and its theoretical errors are discussed. Finally, in section 4 we 
give a brief summary together with some comments on the experimental possibilites. 



2 Hadron Energy Spectrum 

In this section we present the calculation of the hadron energy spectrum which is observed 
in the rest frame of the T(4S') resonance. This frame coincides with the laboratory frame in 
the case of a symmetric 5-factory. We take into account perturbative QCD corrections to 
0{as), but resum Sudakov-like double logarithms. Bound-state effects are evaluated within the 



ACCMM model and the HQET ^ approaches. We anticipate that both methods will lead 
to similar results which, for the present case of a heavy to light transition, depend effectively 
only on one parameter, viz. 6-quark mass. Theoretical accounts as to why the two distinct 
descriptions of bound-state effects in heavy to light transitions should agree, have been given 
in Ref. fTB . 



2.1 Hadron Energy Spectrum in the ACCMM Approach 

In the ACCMM model the 5-meson consists of a 6-quark and a spectator antiquark q, fly- 
ing back-to-back in the B-rest frame with 3-momentum vectors p and —p, respectively; the 
momentum distribution of the spectator is assumed to be of Gaussian form 

= ^ exp ; p=\p^ , (1) 
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normalized according to 

/•oo 

dpp^^{j)) = l . (2) 



The basic feature of this model is the requirement that in the B rest frame the energies of the 
two constituents have to add up to the 5-meson mass m^. This is possible only when at least 
one of the masses of the constituents is allowed to be 3- momentum dependent. Usually, the 
mass of the spectator m^p is considered to be fixed to a constituent quark mass value, while 
the 6-quark mass becomes momentum dependent. 



rribipf = m| + m^^ - 2mB \lv^ + "^sp • (3) 

We now consider the semileptonic decay B = {b,q) u{+g)ii7 + q. The symbols pb, Pu, Pg, 
Pi, Ptj, Pg, Pb, and pn denote the four- momenta of the 6-quark, the w-quark, the gluon, the 
charged lepton, the anti-neutrino, the spectator, the i?-meson, and the hadronic matter of the 
final state, respectively. In addition, we define the four vector g'^ 

q = Pi+Pi7 ■ (4) 

Due to the fact that p^ + Pg = Pb ^^e ACCMM model, the vector q^^ is the same with or 
without the spectator four- momentum, viz., 

q = Pb - PH = Pb + Pq - Pu{-Pg) - Pg = Pb - Pu{-Pg) ■ (5) 

For this reason, it is technically easier to concentrate first on the double differential distribution 
dPV / dq^dq^ , because the spectator does not appear explicitly. To achieve this, we first derive 
the double differential decay rate (including the 0{as) radiative corrections) for h — >• u{+g)tu, 
where the 6-quark decays at rest. For simplicity, the w-quark is assumed to be massless. The 
kinematical variables q^ and g° vary in the region 

< < , 

For each q^, the tree- level process (as well as the virtual gluon corrections) are concentrated at 
the upper endpoint of the g° range in Eq.(|^). In the following it is more convenient to introduce 
dimensionless variables x and y which vary independently in the range [0, 1]: 

2 2 2 

q = X m^, 

gO= [x + i(l-a;)2y]mb. (7) 

With the above change of variables, we now proceed to calculate the double differential 
distribution d'^T/dxdy. In terms of the new variables {x,y), the tree-level contribution and the 
virtual corrections are concentrated at y = 1, while the bremsstrahlung corrections give rise to 
a continuous distribution in the whole {x, y) domain. Individually, both the virtual corrections 
and bremsstrahlung contributions suffer from infrared and collinear divergences, which occur 
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at y = 1 for a given value of x. However, if the variable y is integrated over a range Sq < y < 1 
{sq < 1), these singularities cancel, such that the quantity 



dT 

ax 



1 £T 



dy 



1 d^T 



dy 



-0 d^T 



brems 



dy 



(8) 



/so dxdy Jo dxdy Jo dxdy 
remains finite. The first term on the right hand side (RHS) in Eq.(H) has been calculated in 



Ref. |T^; the result is 

dV fi d^r 



dy = Ax{l-xy{l + 2x') [1 



2a. 



G{x)]T, 



dx Jo dxdy ^ / v ■ ' i 

where Fq and the function G{x) which containes the radiative corrections are given by 
Glml\V^,\' 



(9) 



To 
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(10) 



x] 



[8x2(1 - x2 - 2X^) logX + 2(1 - x2)2(5 + 4x2) log(l - X2) - (1 - x2)(5 + 9x2 _ g^4)] 



2(l-x2)2(l + 2x2) 



+7r2 + 2Li(x2) - 2Li(l - x^) . 
Here Li(x) is the Spence function defined as 



df 

Li(x) = -/ ^log(l-t) 
Jo t 



(12) 



Because the tree-level contribution and virtual corrections are concentrated aty = l, the second 
term on the RHS of Eq. (H) contains the gluon bremsstrahlung contribution only (as indicated 
by the notation). As the endpoint region is cut off by Sq < 1, this term is finite; consequently, 
the infrared and coUinear regularization is not necessary from the very beginning. We have 
calculated this term; the result is 



- fr^^^ ^ ^^^^ _ ^ 2a;2) ^ (log2(l - So) + H{so)) To 

dxdy Sir ^ ' 



(13) 



where 



Jo ^1 — " 



y - 2 

1^1 — x) (3 + X + xy — y) 



log(l -y) -2 log 



2 - ?/(l - x) + K 



:i + x)2 

K r7(l + x)(l + 2x2) 



2(l + x)2(l + 2x2) 



1-y 



;i-x)(3-2x2 



(14) 



The quantity k in Eq. (0) is defined as k = \jy'^{\ — x)^ + 4x?/. Combining the two contribu- 
tions, Eqs.(p|) and (|13|), we obtain 



-1^(.0)=4X(1-X2)2(1 + 2X2) 

1 «x 



2a, 
37r 



log2(l-so) 



2«. 
37r 



(G(x) + ff(so)) 



e(l-so). (15) 



The double logarithms arise from the soft and coUinear gluons and become important as Sq ^ 1. 
Resumming these double logarithmic terms to all orders, we get 



-^^(so) = 4x(l-x2)2(l + 2x2) exp (-^ 
i. ax \ ovr 



log'(l 



So, 



1 - 



2a, 



(G(x) + if(so)) 



e(l-so). 
(16) 

This expression enables us to reproduce the Sudakov exponentiated double-differential decay 
rate by differentiating with respect to Sq. 



1 d^T 
Tq dxdy 



2a, 



X 



d / 1 dV , 

dSo^LfjdX ^ SQ=y 

-4a;(l - x'^ fil + 2x'^) exp ( - log2(l - y) 

/4a, log(l-|/) r 2a, 
\ 3tc {1 — y) L 3tt 



{G{x)+H{y)) 



2a, dH 
3n dy 



(17) 



To get the parton level hadron energy spectrum for a b quark decaying at rest, we first 
re-express Eg . (plTl) in terms of the variables {q'^,q^) and then integrate over g^; this leads to 



the distribution dT/dq^. As the hadronic energy Ehad is related to g° by Ehad = rni, — 
the spectrum dT/dEhad for a 6-quark decaying at rest is readily obtained. To get rid of the 
6-quark mass dependence mf in the decay width (as seen e.g. in Eq.(0)) and errors thereof, we 
present in the following the differential branching ratio dBR/dEhad that is obtained by dividing 
dT/dEhad by the theoretical semileptonic 6-quark decay width and multiplying this ratio 
by the experimentally measured semileptonic branching ratio BRsi = (10.4 ± 0.4)% |[T8[| : 



dBR{B XJ 



dT 



dE. 



had 



r., dE, 



BR 



si- 



had . 



(18) 



The semileptonic decay width, neglecting the small B — > X„£z/ contribution, reads 



si 
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g{mc/mb) 1 - 



2as{mb] 
37r 



f{mc/mb) 



(19) 



where the phase space function g{u) is defined as 



9[u) 



l-8u^ + 8u^ - 



24m logM 



(20) 



and the radiative correction function in an approximate analytic form is given by |TT|: 



fin) 



31 



'l-u) 



(21) 



The result (based on the the Sudakov improved differential branching ratio as just derived) 
is shown by the dash-dotted line in Fig. 1 (using rrih = 4.85 GeV and = 1.45 GeV). To 
emphasize the effects of the Sudakov resummation, we have also plotted the result obtained 
in pure 0{as) perturbation theory |19|, i.e., without exponentiation (short-dashed line). This 
curve exhibits a (integrable) double logarithmic divergence when E^ad = mbf^ is approached 
from above. The effects of the exponentiation are therefore most prominent in the region around 
mft/2 as illustrated in Fig. 1, because Sudakov-exponentiation suppresses the singularity just 
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mentioned. We stress, however, that the effect of exponentiation is neghgibly small in the 
whole region below the charmed hadron threshold we are mostly interested in. This should be 
contrasted with the case of the energy spectrum of the charged lepton near the kinematical 
endpoint where these effect are most pronounced. 

Another interesting aspect concerning the radiative corrections to the hadron energy spec- 
trum is that the kinematical boundary of the hadron energy depends on whether the final state 
contains bremsstrahlung gluons or not. In the absence of bremsstrahlung gluons, the kinemat- 
ical endpoint is at mb/2 = 2.425 GeV as can be seen form the long-dahed curve representing 
the result without QCD corrections. Fig. 1 clearly illustrates that, above E^ad = 2.425 GeV, 
a significant tail due to gluon bremsstrahlung, hence of order as{mb) in strength, is present in 
the spectrum. 




Figure 1: The long-dashed line shows the hadron energy spectrum for a 6-quark decaying at 
rest without taking into account any QCD corrections. The short-dashed line is the corre- 
sponding spectrum including 0{as) virtual and bremsstrahlung corrections. The result after 
exponentiating the Sudakov double logarithms (discussed in the present section) is shown by 
the dash-dotted line. The solid line shows the hadron energy spectrum for a 5-meson decaying 
at rest. The bound state effects are calculated with the ACCMM model with pp = 344 MeV 
and rrisp = 150 MeV (corresponding to (rrib) = 4.85 GeV). The open circles show the hadron 
energy spectrum due to a 5-meson decaying at flight with a momentum |pb| = 330 MeV. 

We now turn to implement the bound-state effects using the ACCMM model. We start 
from the double differential distribution dV/dq^dq^ in Eq. ([T^. As the 6-quark moves with 
momentum p inside the S-meson, we first replace the mass nib by mb{p) as given in Eq. (^. 
We then Lorentz boost the double differential distribution and get the spectrum for a 6-quark 
decaying at flight (momentum p). Finally we convolute the spectrum with the ACCMM distri- 
bution function given in Eq. (j^). This leads to the ACCMM averaged quantity dT/dq^dq^ for a 
5-meson decaying at rest. Again, this distribution is straightforwardly converted to dT/dEhad- 
To get the differential branching ratio, we replace the 6-quark mass rrib in the expression for the 
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semileptonic decay width in Eq. ([T9|) by the ACCMM average value {rrik) derived by applying 
the ACCMM convolution to the total semileptonic decay rate. This ACCMM averaged 6-quark 
mass is given by 

{ml) = I dpp^{mbip)f^p) ■ (22) 

The rrih — rric mass difference is reliably calculated to be 3.40 GeV by HQET pO|. Therefore, 
we set rric = (mb) — 3.40 GeV in the Tgi- The result is shown in Fig. 1 by the solid line; here we 
have used pp = 344 MeV and m^p = 150 MeV, which yields (rrib) = 4.85 GeV according to Eq. 
(|22|) . As shown in Fig. 1, the main effect of the Fermi motion of the 6-quark inside i?-meson 
is to shift the perturbative hadron energy spectrum uniformly to higher energies by about 300 
MeV. We will elaborate this point in more detail when discussing the HQET approach to the 
bound-state effects, where the same result can be understood in a more transparent way. 

There is one more source of Doppler shift to the hadron energy spectrum. At a symmetric 
5-factory the S-meson is produced from the decay of the T(4S') resonance. Due to the released 
binding energy the 5-mesons are produced with a momentum \pb\ ~ 330 MeV in the T(4S') 
rest frame. We have also worked out this effect to the spectrum; the corresponding result is 
shown in Fig. 1 by open circles. As this effect is very small, we will not consider it any more 
in foregoing discussions. 

Although the ACCMM model has two input parameters, viz. pp and rrisp, it turns out that 
the hadron energy spectrum is sensitive only to one parameter, namely, the average 6-quark 
mass {rrib), which is a function of pp and m^p (see Eq. (|22|)) To illustrate this point, we 
have chosen two different pairs of {ppi^sp) values, which both correspond to the same value 
of {rrib) (=4.85 GeV in the present case). The correspondig result plotted in Fig. 2 indicates 
little dependence on pp and m^p, once {mb) is hold fixed. In Fig. 3, we have varied {mb) over 
the range indicated in the plot. We have found that the hadron energy spectrum depends on 
{mb) rather strongly, especially in the region below the charmed hadron threshold. 

2.2 Hadron Energy Spectrum in the HQET Approach 

Next we calculate the hadron energy spectrum using the HQET approach. Since the derivation 



is discussed in detail in [15|, we only repeat the results relevant for our discussions. In the 



following we take into account the leading corrections in and 1/mb, but we neglect mixed 
higher-order corrections of the order a^ — . We assign the b quark the same four velocity as 



the heavy meson and define the dimensionless parton level quantity Eq as 



V ■ {Pb- q) _ ^_ q ■ V 



Eo= "^" =1--—^ . (23) 
mb mb 

The physical hadron energy Ehad is given as 

Ehad = V ■ (Pb - q) = me - q ■ V . (24) 

The leading order equality between msE^ and Ehad gets modified by corrections hnear in 1/mb 

T A1 + 3A2 / ^ Ai + 3A2\ - , , 

Enad = A - ^- ^ + - A + ^0 + ■ ■ ■ , 25 

ImB \ 2mB J 



^ It has been pointed out that transitions from heavy quarks to hght quarks with masses niq < (Amb)^/^ 
always give rise effectively to a one-parameter dependence p6l. 
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Figure 2: The hadron energy spectrum (based on the ACCMM model) is shown for two different 
pairs of model paramters {pp, nisp). The solid line corresponds to (374 MeV, 0) while the dashed 
line represents the result for (252 MeV, 300 MeV). Both set correspond to the same value of 
the 5-quark mass {rrih) = 4.85 GeV. 



where we have used the relationship 



J Ai + 3A2 

niB = mb + A h 

Zrub 



(26) 



between the b quark and the 5-meson masses. The ellipses in Eqs. ( pSf ) and (^) stand for 
higher order corrections in the l/mt expansion. We split the differential spectrum into two 

dEo 

parts: 



dT' 



+ . (27) 

dEo dEo dEo 

The first term on the right hand side (RHS) of Eq. (^) consists of the 1/ml corrections and 
has been calculated in For m„ = one gets (0 < Eo < 1/2) 



1 rfr^/^ 



To dEo 



16E2[3-4i^o] 



+ 16 



Ai 
2m| 



6^0 + 12E^ + — 43 j + (^_3 _ + 36^2 + 2043 



A 



+ 



' 5Ai 
3ml 



m 



S{Eo - 1/2) 



B , 



Ai 



5\Eo-l/2) + 0{as;llml 



(28) 



with To given in Eq. (jT^). The second term on the RHS of Eq. ( P?] ) contains the perturbative 
0{as) corrections; it has been calculated by Czarnecki et al. [|19|. For completeness, we list 
their result for a massless -u-quark: 

1 dV^' ?,2n 

■'g,{Eo) , (29) 



To dEo 



37r 
8 




Figure 3: The hadron energy spectrum based on the ACCMM approach is shown for different 
values of the 6-quark mass (mt) . 



El 


(1 

190 


+0 





TT 



— + Li(l - 2Eo) 



where the function Qi for Q < Eq < 1/2 reads 
Gi{E^) = 

+ (8^0 - 9) log 2^0 + 2(4^0 - 3) 
while for 1/2 < < 1 

^i(^o) = (1 - ^o) (32^0^ - 136io' + 1034^0^ - 2946^^ + 1899^0 + 312) 
180 ^ ' 

log(2Eo - 1) [meI - A^El - 2AEo - 5 



+El (3 - AEo) 



TC 



4Li(l/(2Eo)) + log'(2Eo - 1) - 21og'(2io 



(30) 



(31) 



1 dr 



Using relation (^5|) between Eo and Ehad, the corresponding hadron energy spectrum ^ 



is readily obtained. While the quantity ^^-j§r is free from corrections linear in l/rrib, there are 



such corrections in the hadron energy spectrum The differential branching ratio in the 

HQET approach is then given by 



dBR{B ^ X^eu) _ \Vub\' (i^ 



dE, 



had 



I Kb 1 2 g{mc/mb) 



BR. 



1 z f{mc/mb) + 



1 -1 



37r 



2ml 



, (32) 



where the phase space function g{u) and the radiative correction function f{u) are listed in 
Eqs. ( pOD and (pi]), respectively. The function h{u) contains the l/ml corrections to the semilep- 



tonic decay B — ^ Xciiy; it reads 

h(u) = Ai + \-9 + 24u^ - 72m^ + 72u^ - 15m® - 72u^ log u] . (33) 
g{u) L J 

To illustrate the qualitative features of the non-perturbative corrections, it is instructive 
to retain only the dominant terms which are linear in l/mi,. The relation between Ehad and 



Eq (given in Eq. (p5D) then simplifies to Ehad = A + rribEo = A + Eq, where Eq is the hadron 
energy at the parton level. Therefore, the dominant bound-state effect is simply a uniform 
shift of the parton level hadron energy spectrum to higher energies by an amount of A; exactly 
the same feature was also noted when implementing the bound state effects using the ACCMM 
approach (see Fig. 1). From these considerations it is clear that the most important parameter 
is A ^ rriB — mb, while the parameters Ai and A2 play a relatively minor role. In the numerical 
evaluations we therefore have used the fixed values Ai = — 0.5 GeV^ and A2 = 0.12 GeV^ and 
presented the results as a function of rrib (or, equivalently. A) only. 

In Fig. 4 we compare the ACCMM and the HQET approaches. As an example we have 
chosen rub = (mb) = 4.85 GeV. The plot illustrates that the two methods agree with each other 
rather well, especially in the energy range 1 GeV < Ehad < ttld. We have checked that this 
agreement holds everywhere in the relevant b quark mass range nib = (4.8 ± 0.2) GeV. 




Figure 4: The hadron energy spectrum is shown for two different descriptions of the bound 
state effects. The solid curve shows the result obtained by using the ACCMM model while the 
dashed line shows the HQET result. 



3 Extraction of Vub from the Hadron Energy Spectrum 

As discussed in the Introduction, theoretical predictions for the hadron energy spectrum are 
expected to be rather reliable in the relevant kinematical region below the D-meson mass, once 
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a sufficiently high lower-cut in the hadron energy is made at the same time. Choosing the lower 
cut at 1 GeV, the invariant hadronic mass ranges from up to E^ad > 1 GeV for all values 
of the hadronic energy in the window 1 GeV < E^ad < '>tid ■ In this case a large number of 
different hadronic final states contributes to the spectrum; hence, quark-hadron duality should 
hold quite well. We therefore propose to measure experimentally the kinematical branching 
ratio defined by 



BR^-{B XM = r dBR{B -^X^l.^ ^^^^^ 

J I GeV dEhad 

A comparison with the corresponding theoretical quantity then allows to extract the CKM 
ratio iKtbl/lKftl- While leaving aside the experimental question how accurately the kinematical 
branching ratio can be measured now or in the future, we would like to point out in the following 
that the theoretical uncertainties in this observable are small enough to reduce substantially the 
present theoretical error on \Vub\- As it turns out that the theoretical error of the kinematical 
branching ratio is dominated by the uncertainty of the effective b quark mass (i.e., in the 
HQET approach and {rrih) in the ACCMM approach), we take into account only this effect 
in the following analysis 0. A reasonable range for the value of rrif, can be inferred from the 
measurement of the lepton energy spectrum in the inclusive decays B Xc^^iu. Fitting to the 
ACCMM model, CLEO extracted pp to be 264 ± 16 MeV (using a value of rusp = 150 MeV 
for the mass of the spectator quark) p|, |2ll]. In a more recent analysis using lepton tags |]22 



they extracted the value Pf = 347 ± 68 MeV. We thus choose to vary pp m the somewhat 
larger range 200 MeV < Pf 435 MeV which covers both measurements. Using Eq. (p2D , this 
Pi?-range translates into the (m;,)-range 

4.75 GeV < {nib) < 5 GeV. (35) 

We note that this range is also compatible with the best value ((m{,)=4.77 GeV) fitted from 
the photon energy spectrum in B ^ Xg'y . 



So far, we have not discussed the question, if the cascade decay b — > cX siuX , where 
the symbol X denotes light quarks, can fake a 6 — uiu transition. If the c-quark is off-shell, 
this process can kinematically fake a 6 — uiu transition, but in this case the process is higher 
order in the weak interaction, hence negligible. On the other hand, if the c quark is on-shell, 
the energy of the hadronic system X is larger than m^). Therefore, this cascade process is not 
a backgound to the b —>■ uiu transition; this is in contrast to the endpoint analysis of the lepton 
energy spectrum, where the cascade process has to be subtracted out. 

To get an idea what fraction R of the total semileptonic b ^ u events will be captured in the 
energy window 1 GeV < Ehad < ""^d, we have used the present central value for |\46|/|Kb| = 0.08 
[p^ . The mfe dependence of the kinematical branching ratio and of the fraction R are shown in 
Table 1 for both the ACCMM and the HQET approaches. Varying rrib in the range specified in 
Eq. (^3]) and using the results for BR^q^j^ in Table 1 (where the variation is somewhat larger 
than in the ACCMM model), we obtain 



cb 



^The dependence of the kinematical branching ratio on the numerical value of is very small. In our 
numerical evaluations we have taken as — 0.205. 
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rrih (GeV) 


B^ACCMM * 10'' 


Raccmm (%) 


BRTqet * 


Rhqet{%) 


4.60 


2.41 


23 


2.17 


21 


4.65 


2.58 


24 


2.40 


22 


4.70 


2.76 


25 


2.64 


24 


4.75 


2.94 


26 


2.89 


25 


4.80 


3.15 


27 


3.14 


27 


4.85 


3.37 


28 


3.41 


28 


4.90 


3.59 


29 


3.68 


29 


4.95 


3.83 


30 


3.95 


31 


5.00 


4.07 


31 


4.23 


32 



Table 1: The kinematical branching ratio BR^^^ defined in Eq. (0) and the corresponding 
fraction R of the semileptonic h ^ u events lying in the energy window 1 GeV < E^ad < f^D 
are given as a function of m;, for the ACCMM and the HQET approach. |Kjb|/|V^b| = 0.08 is 
assumed. 



Denoting the measured kinematical branching ratio by BR^!1, one can extract |\4fe|/|V^f,| to be 



"exp 



Wr 



cb\ 



BRl^ ■ (4.30 ± 0.41) . (37) 



This implies that the theoretical error of the ratio is approximately ±10%. Taking 

the somewhat larger range m;, = (4.80 ± 0.15) GeV adopted in [|T], we get from Table 1 



IK 



ub 



Wr 



cb\ 



BR^il ■ (4.60 ± 0.56) , (38) 



which implies a theoretical error of about 12 %. 

To illustrate that our proposal has the potential to lead to a more precise determination 
of Vub, it is instructive to briefiy review the present situation. The traditional inclusive lepton 
endpoint spectrum analysis done at CLEO and at ARGUS leads to the result = 
0.08 ± 0.03, where the error is dominated by theory |2^. A recent new input to this quantity 
is provided by the measurements of the exclusive semileptonic decays B — >• (tt, p, a;)£z/ ||. 
The value for \Vub\ quoted in the most recent analysis |^ is \Vub\ = (3.3 ± 0.2to;4 ± 0.7) x 10^^, 
where the errors are statistical, systematic and estimated model dependence, after excluding 
models which are unable to predict the correct ir/p ratio. 



4 Summary and Discussion 

In this paper we have proposed to measure the hadron energy spectrum from inclusive semilep- 
tonic B decays in order to extract the CKM matrix element Vub with improved precision. The 
main advantage of our proposal is that the energy window, in which the b —>■ c£~V transition 
is kinematically absent {Ehad < ^d), is much wider than the corresponing window available 
in the traditional lepton energy endpoint spectrum analysis. Even after imposing a relatively 
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high lower-cut at Ehad = 1 GeV (in order to avoid the region of phase space where the range 
in the invariant hadronic mass is too narrow to invoke quark-hadron duahty), a much larger 
fraction of the b — >■ ui~V events is captured in the remaining window 1 GeV < Ehad < "mn than 
in the lepton spectrum endpoint analysis. 

After calculating the hadron energy spectrum in QCD improved perturbation theory, we 
have implemented the bound-state effects using both the ACCMM and the HQET approaches; 
the two methods gave essentially the same result, in particular in the relevant kinematical 
window. Qualitatively, the dominant bound state effect is a uniform shift of the parton level 
spectrum. We have pointed out that the theoretical error in the ratio iKibl/lKfel; which is 
dominated by the present uncertainties in the rrih quark mass, is about a factor of 2 to 3 
smaller than the one in the lepton endpoint spectrum analysis and about a factor of 2 smaller 
than the model uncertainties in the branching fractions of the exclusive decays B —>■ (vr, p, uj) iu. 
Therefore, for both statistical and theoretical reasons, the extraction of \Vuh\ from the hadron 
energy spectrum seems very attractive. 

On the experimental side, our proposal is most suited for a symmetric i?-factory running 
at the T(4S') resonance, which is currently available at the CLEO experiment. Tagging the 
events from the T(45') decay, in which one i?- meson is decaying semileptonically and the other 
one nonleptonically, the energy of the final state hadrons stemming from the semileptonically 
decaying 5-meson is easily obtained by adding up the energies of all the hadrons in the final 
state and then subtracting mr/2. On the other hand, in asymmetric B factories, the hadron 
energy spectrum is harder to measure because one first has to reconstruct the corresponding 
distribution in the rest frame of T(45'); to perform the corresponding boost, one has to measure 
precisely both the energy and the momentum of each final state hadron, which requires accurate 
particle identification. 

The spectrum of the invariant hadronic mass {dV/dmx) in inclusive semileptonic decays 
B Xc^u^jy is another source from which one may try to extract Vub- Requiring mx to be 
below m/), the process B Xciiy can be totally supressed. Consequently, the integral of the 
hadron invariant mass distribution below is another observable proportional to iKitp. As 
the invariant mass is integrated over a large range, which covers about 95 % of all the b — *• ui~u 
events [^, this observable is also theoretically viable. However, we believe that the invariant 



mass spectrum is more difficult to measure than the hadron energy spectrum proposed in this 
paper, because first, one has to measure the four momenta of all the final state hadrons from 
the T(4S') decay and second, one has to find a subset of final state hadrons with an invariant 
mass of m^, corresponing to the 5-meson which decays hadronically. Only then the invariant 
mass of the semileptonically decaying 5-meson can be determined. 

Finally, we point out a potentially interesting possibility of a direct measurement of as{mb) 
or A. We have shown that, once the real gluon bremsstrahlung correction is taken into account, 
the kinematic maximum of the hadron energy shifts from the tree-level (and virtual gluon 
correction) endpoint (m^ + mg)/2mb to mb- For the b ^ u case. Fig. 1 shows that the 
bremsstrahlung tail extends approximately 400 MeV beyond mb/l. We also have shown that 
the dominant bound-state effect is to shift the spectrum by A for both ACCMM and HQET 
approaches. Therefore, once A is known accurately, then as{mb) can be extracted directly from 
the bremsstrahlung tail spectrum. Alternatively, given an accurate value of ^^(^b); the non- 
perturbative parameter A can be accurately measured. Since the bremsstrahlung spectrum of 
b ^ c extends further out than that of 6 — > m, the above proposal should be most suited for the 
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b ^ c transition. In this case, one has to take into account the finite c-quark mass dependence 
of the perturbative QCD corrections and bound-state effects. 

We thank A. Ah, S.J. Brodsky, B. Grinstein, S.-K. Kim, J.H. Kiihn, Y.-J. Kwon, M. Lu, P. 
Minkowski, R. Pohng, E. Thorndike, and R. Wang for helpful discussions. 
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